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We shall investigate the asymptotic properties of the Bianchi type IX cosmological model in the 
Brane-world scenario. The matter content is assumed to be a combination of a perfect fluid and 

> 

t^J- ■ a minimimally coupled scalar field that is restricted to the Brane. A detailed qualitative analysis 

ON ' 

of the Bianchi type IX braneworld containing a scalar field having an exponential potential is 

(N 

undertaken. It is found that the Brane- Robertson Walker solution is a local source for the expanding 

(N ' 

Bianchi IX models, and if k 2 < 2 the expanding Bianchi IX models asymptote to the power-law 
inflationary solution. The only other local sink is the contracting Brane-Robertson Walker solution. 

Sh '. 

(3JT), An analysis of the Bianchi IX models with a scalar field with a general potential is discussed, and 

it is shown that in the case of expanding models, for physical scalar field potentials close to the 
initial singularity, the scalar field is effectively massless, and the solution is approximated by the 
Brane-Robertson Walker model. 
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I. INTRODUCTION 



It is believed that Einstein's General Relativity breaks down at sufficiently high energies. 
Developments in string theory suggest that gravity may truly be a higher dimensional the- 
ory, becoming an effective 4-dimensional theory at lower energies. Some researchers have 
suggested an alternative scenario in which the matter fields are restricted to a 3-dimensional 
Brane- world embedded in 1 + 3 + d dimensions (the bulk), while the gravitational field is 
free to propagate in the d extra dimensions |]T]. In this paradigm it is not necessary for 
the d extra dimensions to be small, or even compact, a radical departure from the standard 
Kaluza-Klein scenario. Randall and Sundrum have shown that for d = 1, gravity can be 
localized on a single 3-brane even when the fifth dimension is infinite. It has now become 
very important to test the astrophysical and cosmological implications of these higher di- 
mensional theories derived from string theory. Can these cosmological models derived from 
string theory explain the high degree of homogeneity and isotropy we currently observe? 

The dynamical equations on the 3-brane || |], || differ from the general relativity equa- 
tions by terms that carry the effects of imbedding and of the free gravitational field in the 
five-dimensional bulk. The local (quadratic) energy-momentum corrections are significant 
only at very high energies and the dynamical equations reduce to the regular Einstein field 
equations of General Relativity for late times. However, for very high energies (i.e, early 
times), these additional energy momentum correction terms will play a very critical role 
in the evolutionary dynamics of these Brane-world models. In addition to the matter field 
corrections, there are nonlocal effects (that modify the Friedmann equation on the brane in 
these models) from the free gravitational field in the bulk, transmitted via the projection 
Efjtu of the bulk Weyl tensor, that contribute further corrections to the Einstein equations 
on the brane. 

Cosmological observations indicate that we live in a universe which is remarkably uniform 
on very large scales. However, the spatial homogeneity and isotropy of the universe is difficult 
to explain within the standard general relativistic framework since, in the presence of matter, 
the class of solutions to the Einstein equations which evolve towards a FRW universe is 
essentially a set of measure zero ||. In the inflationary scenario, we live in an isotropic 
region of a potentially highly irregular universe as the result of an expansion phase in the 
early universe thereby solving many of the problems of cosmology. Thus this scenario can 
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successfully generate a homogeneous and isotropic FRW-like universe from initial conditions 
which, in the absence of inflation, would have resulted in a universe far removed from the one 
we live in today. However, still only a restricted set of initial data will lead to smooth enough 
conditions for the onset of inflation (i.e., the so-called cosmic no-hair theorems only apply 
to non-generic models - see 0, so the issue of homogenization and isotropization is still 
not satisfactorily solved. Indeed, the initial conditions problem, that is to explain why the 
universe is so isotropic and spatial homogeneous from generic initial conditions, is perhaps 
one of the central problems of modern theoretical cosmology. These issues were recently 



revisited in the context of brane cosmology |13|, [L4J], and this is one of the motivations for 
the present work. 

Indeed, researchers have investigated both anisotropic and isotropic Brane-world models, 
trying to ascertain whether the effects of the bulk gravitational field would allow one to 
solve the isotropy problem. A lot of effort has been directed at the so-called Friedmann 
Brane-world models || [T(J. Dynamics of a brane-world universe filled with a perfect fluid 
have been intensively investigated during the last three years fll , 12 , |13| , |14|| . It has been 
found there exists new regimes that are not inherent in the standard cosmology, such as 
stable oscillation ]12| and the collapse of a flat universe [15|]. Some features of brane-world 
inflation have been studied in ||, W, T7] and the cosmo logical dynamics for exponential 
scalar field potentials have been described in [ |T8| , |H) |. 

Shiromizo et al. || have developed an elegant covariant approach to the bulk effects 
on the brane. The equations derived by Maartens Q are an extension of earlier work by 
Ellis and MacCallum that has been subsequently developed more recently in the book by 



Wainwright and Ellis p0[ . Using the formalism developed by Maartens and Wainwright 
and Ellis, we propose to investigate the dynamical behavior in a wider class of anisotropic 
models than what has been previously analyzed. Spatially homogeneous models of Class A 
(and in particular Bianchi type IX), containing both a perfect fluid and a scalar field will 
be investigated. The resulting field equations will yield a system of ordinary differential 
equations, suitable for a geometric analysis using dynamical systems techniques. This anal- 
ysis will determine whether the dynamics of the Brane-world scenario mimics the dynamics 
of a General Relativistic cosmology at late times. We are particularly interested in both 
the early-time (nature of the initial singularity) and late-time behaviour (i.e., whether these 
models inflate and isotropize). 
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II. GOVERNING EQUATIONS 



The field equations induced on the brane are derived via an elegant geometric approach 
by Shiromizu et al. 0, ||, using the Gauss-Codazzi equations, matching conditions and Z 2 
symmetry. The result is a modification of the standard Einstein equations, with the new 
terms carrying bulk effects onto the brane: 

G = —Ag^y + k T^ u + k S^ u — E^y , (2-1) 

where 



2 _ 8?r \ _ q ^ ^ 4:71 



A, |^| 

3M p 3 



(2.2) 



It is common to assume through fine tuning (a la Randall Sundrum) that the effective 
cosmological constant on the brane is zero, i.e., A = 0, However, we shall assume that it is 
non-zero but positive. 

The brane energy-momentum tensor for a perfect fluid and a minimally-coupled scalar 
field is given by 

rp rp perfect fluid < rp scalar field (n o\ 

where 

^perfect fluid = + ^ ^ 

^scalar field = ^ _ g ^ + V ^)\ , ( 2 . 5 ) 

where is the fluid 4-velocity, p and p are the energy density and isotropic pressure, 
is the minimally coupled scalar field having potential V ((/>), and the projection tensor 
h^ v = g^y + u^Uy projects orthogonal to u^. If </> ;M is timelike, then a scalar field with 
potential V{<p) is equivalent to a perfect fluid having an energy density and pressure 

p scalar field = (f,.^ + ( 2 .6) 

p scalar field = __^JM _ y(0). (2.7) 

The bulk corrections to the Einstein equations on the brane are of two forms: firstly, 
the matter fields contribute local quadratic energy-momentum corrections via the tensor 
S^y, and secondly, there are nonlocal effects from the free gravitational field in the bulk, 



4 



transmitted via the projection of the bulk Weyl tensor, 8^ v . The local matter corrections 
are given by 

Sfj,v = —T a a T^ u — -T^T 01 ^ + —g^iv [3T a/3 T al3 — (T a a ) 2 ] , (2.8) 



which is equivalent to 



g perfect fluid = + L p {p + 2 p) h, u , (2.9) 



for a perfect fluid and 



i(-i, 

6 V 2 

for a minimally-coupled scalar field. If we have both a perfect fluid and a scalar field and if 
we assume that the gradient of the scalar field <p' ,fl , is aligned with the fluid 4- velocity, u M , 



that is \J —4>;a4> ,a = u>1 - ( m general ;M need not be aligned with thereby creating a 
rich variety of cross terms.) The local brane effects due to a combination of a perfect fluid 
and a scalar field are then 

+ I2 ( P " + V )[p + 2 P- \<f>;«<P a ~ W)) V (2- 11 ) 

The non-local effects from the free gravitational field in the bulk are characterized by 
the projection of the bulk Weyl tensor onto the brane. Given a timelike congruence on the 
brane, the bulk correction, £^ v can be decomposed [|J via 

4 r i i 

(2.12) 



'-'[11' 



U{upU v + ^V-) + TV + 2 Q(v u v) 



(See [|J for further details.) In general, the conservation equations (the contracted Bianchi 
identities on the brane) do not determine all of the independent components of £ M „ on the 
brane. In particular, there is no evolution equation for TV and hence, in general, the 
projection of the 5-dimensional field equations onto the brane does not lead to a closed 
system. However, in the cosmological context studied here, we will assume 

D M W = g M = 7V = 0, (2.13) 
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where D M is the totally projected part of the brane covariant derivative. Since = 0, in 
this case the evolution of E^ v is fully determined |2j . In general U = U(t) ^ (and can be 
negative) in the Friedmann background || [Hj. For a spatially homogeneous and isotropic 
model on the brane, equation ( |2.13| ) follows, and similar conditions apply self-consistently 
in other Bianchi models ||. In the appendix we prove that the integrability conditions for 
V^y = 0, = implies spatial homogeneity. Physically V^v corresponds to gravitational 



waves and will not affect the dynamics close to the singularity ||25|| . From the analysis of 
the evolution equation for close to the initial singularity, it can be shown that a small 
does not affect the dynamical evolution of U (to lowest order) f|, [13|, [26| . 
All of the bulk corrections mentioned above may be consolidated into an effective total 
energy density and pressure as follows. The modified Einstein equations take the standard 
Einstein form with a redefined energy-momentum tensor: 



G 



fJ,U 



2rT-itotal 
K 1 [IV 



where 



rptot — 



A it u*l<7 Lf 



(2.14) 



(2.15) 



is the redefined equivalent perfect fluid energy-momentum tensor with the total equivalent 
energy density due to both a perfect fluid and a scalar field 



P 



total 



A 



K- 



+ P + 



£ 4 
K 



(2.16) 



total 



P 



A 

K 



+ P + 



V(0) 



^(p-\^ a + V{<P) \ (p + 2p-| 



l <P a -v(4>) )+hi 



(2.17) 



where we have assumed that D^U 



Qu = v, 



in the cosmological case of interest here. 



As a consequence of the form of the bulk energy-momentum tensor and of Z2 symmetry, 
it follows 0] that the brane energy-momentum tensor separately satisfies the conservation 
equations, (where we have tacitly assumed that the scalar field and the matter are non- 
interacting) i.e., 

^perfect fluid = q which yieldg ^ + 3^ + p) = (2.18) 

• dV 

which yields <j> + 3H<p + — = , (2.19) 



rp/j, scalar field 



whence the Bianchi identities on the brane imply that the projected Weyl tensor obeys the 
constraint 

S^. fl = K' i S^ ll which yields U + 4HU = (2.20) 

III. BIANCHI IX MODELS 

A. Setting up the Dynamical System 

We shall use the formalism of Hewitt, Uggla, and Wainwright introduced in [^0] for 
positive curvature models (see pages 179-182 in |2(|). The source term (restricted to the 
Brane) is a non-interacting mixture of ordinary matter having energy density p, and a 
minimally coupled scalar field, 0. We shall assume that the matter content is equivalent to 
that of a non-tilting perfect fluid with a linear barotropic equation of state for the fluid, i.e., 
p = (7 — l)p, where the energy conditions impose the restriction p > 0, and the constant 
7 satisfies 7 G [0,2] from causality requirements. We shall also assume that the scalar field 
potential has an exponential form V((f)) = Voe kK ^ = Vqb M l k ^ |22] ; [23], 2l[ . The energy- 



momentum tensor describing the source is given in equation ( |2.3| ) where, for a homogeneous 
scalar field, = (j>(t). 

Our variables are the same as those introduced by Hewitt, Uggla and Wainwright (see 
page 180 in [EIJ), with the addition of 

~ 3D 2> ^ "A - 3D2 , 9- 3D2 , 

where 



D = \lH 2 + J(nm 2 n 3 ) 2 /3. 



The total equivalent dimensionless energy density due to all sources and bulk corrections is 

K "i total _ n 2 /_~_\2 

Q total = ^-^ = Q + Q A + ^ + ^ + a 2 Q u + -D 2 (Q + ^ + ^ 2 j , (3.2) 



where 

2 _ K 



a 

The governing differential equations for the quantities 

X = [D, H, E 1} S 2 , N u N 2 , N 3 , Q, Q A , $, Q u ] 



are as follows 



D' 


= -(l + q)HD 


(3.3a; 


H' 


= q{H 2 - 1) 


(3.3b; 


K 


= (q-2)Ht + -S + 


(3.3c; 




= (q - 2)Ht_ - 


(3.3d; 




= N 1 (Hq-4E+) 


(3.3e; 


K 


= iV 2 (iJg + 2S + + 2v / 3S_) 


(3.3f; 


K 


= N 3 (Hq + 2t + -2V3E_) 


(3.3 g ; 


rV 
s I 


= n\l {2{q + i) — 6^j) 


(^o.ori / 


rV 


= 2f/i2 A (q + 1) 


(3.3i y 


q?> 


= (q-2)H^ -^yk$ 


(3.3j; 


$' 


= 2$ (^(1 + q)H + y^-kVj 


(3.3k; 




= 2M u {q-l) 


(3.3i; 



The quantity q is the deceleration parameter, and S + and SL are curvature terms that 
are defined by the following expressions: 
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2Y? + 2£ 2 | ^ ^ r > ° 1 ° ,T ' 2 



-fi - fi A + 2^ 2 - $ 



fi + ^ 2 + $) ((37-l)fi + 5§ 2 
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+a 2 fi w + —D 



~ (iV 2 -iY 3 ) 2 -^iVx (2 A, Ao A', 



S_ = l - V3 (N 2 - iV 3 ) (-JVi + N 2 + JV 3 ) 



(3.4) 

(3.5) 
(3.6) 



The evolution equations for p.3h| ) comes from the conservation equation ( j2.18j ). The 
evolution equations for ( |3.3j| ) and ( ^.3k| ) are derived from the Klein-Gordon equation de- 
rived from the conservation equation ( |2.19| ). The evolution equation ( |3.31| ) comes from the 



conservation equation fl2.20|) . 

In addition there are two constraint equations that must also be satisfied, 

Gi(X) = H 2 + ^(iVxiV^) 2 / 3 -1 = 



G 2 (X) 



t% - S 2 
8 



n 



total 



v = o 



(3.7a) 
(3.7b) 



where 



(iV 2 + Nl + Nl - 2{N 1 N 2 + iVxiVg + iy 2 iY 3 ) + ^N^Ns) 2 '* 



Equation Q3.7a ) follows from the definition of D, and equation (|3.7b|) is the generalized 



Friedmann equation. We now have determined the equations describing the evolution of the 
Bianchi type IX brane world models. The resulting equations are suitable for a qualitative 
analysis using techniques from dynamical systems theory. In general, the system of equations 
( p. 3D can be interpreted as X' = F(X) where F : X G IR 12 — > R 12 . We must also make careful 
note of the two constraint equations Gi(X) = and G 2 (X-) = 0. These constraint equations 
essentially restrict the dynamics of the dynamical system X' = F(X) to lower dimensional 
surfaces in R 12 . In principal, these constraint equations may be used to eliminate one of the 
twelve variables provided the constraint is not singular. 

B. Symmetry Transformations 



Since the dynamical system ( |3.3j ) is invariant under the transformation $ — ► — <3> we can 
restrict our state space to D > [by definition of D, see Eq. (|3.1|)] and $ > without loss 
of generality. 

Note the dynamical system (|3.3|) is also invariant under the transformation 

(£+, N u N 2 , N 3 ) -> - ^E_, - N 2 , N 3 , N x ). (3.8) 

This symmetry implies that any equilibrium point with a non-zero S-t term, will have two 
equivalent copies of that point located at positions that are rotated thru an angle of 27r/3 
and centered along a different axis of the N a . 



C. Invariant Sets 

If we assume the weak energy condition for a perfect fluid (i.e., p > 0), then we must 
restrict the state space to the set of points Q > 0. Since we are investigating the behaviour 
of the Bianchi Type IX brane world models in particular, we can restrict the state space to 
N a > without loss of generality. Therefore the state space for the Bianchi IX brane world 
models is the set of points S = {X G R 12 |Gi(X) = 0, G 2 (X) = 0, Q > 0, N a > 0, D > 0, 1 > 
0}. 



The evolution equation for flu, ( P-31Q implies that the surface flu = divides the state 
space into three distinct regions, U + = {X G «S|fi^ > 0}, W° = {X G <S|O w = 0}, and 

w- = {xg «s|ft w < o}. 

There are various invariant sets associated with the matter content. We define six Matter 
invariant sets as 

\fl = 0,$ = 0,tf = 0}, 
\fl = o,$ = 0,§ ^ 0}, 

|n = o,$ ^ o,^ ^ o}, 
|n ^ o,$ = o,^ = o}, 
\n ^ o,$ = o,^ ^ o}, 

|O^0,<M0,^0}, 
where the notation is interpreted as 

(value of ^^(value of i, value of 

In addition to the matter invariant sets, there are invariant sets associated with the 
geometry of the spacetime. The dynamical system ( |3.3| ) implies that any combination of 
the conditions N a > and Np = defines an invariant set. Since these conditions yield the 
Bianchi type of the underlying geometry we call these invariant sets, Bianchi invariant sets. 
The Bianchi type corresponding to the different combinations are 



o o,o 


= {X G S 




= {X G S 




= {XeS 




= {XeS 




= {XeS 


+fl+'± 


= {XeS 



B(I) = {X e S\Ni = 0,N 2 = 0, N 3 = 0} 

B(II) = {X e S\N! ^0,N 2 = 0, N 3 = 0} 

B(VII ) = {X e S\Nt = 0,N 2 ^ 0, N 3 ^ 0} 

B(IX) = {X e Sl^ ^ 0, N 2 ^ 0, N 3 ^ 0} 



where the sets B(II) and B(VIIo) have two additional and equivalent disjoint copies of 
themselves that are determined by the transformation (|3.8|). 



From ( |3.7a| ) we have that — 1 < H < 1 and NiN 2 N 3 < 8. Furthermore, in the invariant 
sets U + and U°, using equation ( |3.7b| ) it can be shown that 



< t 2 + ,t 2 _ ,fl,fl A ,ty 2 ,$,V < 1. 
10 



However, knowing that < V < 1 and < N1N2N3 < 8 is not sufficient to place any bounds 
on the N a 's or D. Furthermore, in the invariant set U~ , we cannot place upper bounds on 



any of the variables without some redefinition of the dimensionless variables ([371 



D. First Integrals 

It is possible to show that the function 

W = (H 2 - 1)[- 2 - 3 ^ ^+^- 2 M fi» % (3.9) 

(where s, t are parameters that can take any value) is a first integral of the of the dynamical 
system ( |3.3|) , that is W = 0. For the particular values s and t we obtain the following 
invariants (for any value K), 

( s = 0, t = 1), k(h 2 - 1) 37 = n ( l y ~ 2) n 2 (3.10) 

( s = 1, t = 0), k(h 2 - 1) 2 = n A n u (3.11) 

( s = -3 7 ,t = 2), KU? U < = fii 37 " 4) fi 4 (3.12) 

(a = _(3 7 _2),t = 1), K^~ 2) = (H 2 - lf^Cl 2 (3.13) 

In the invariant set U + UU°, it is possible to show that q > —1, which implies that — > 
as r — > 00. Using the first and second invariants, we easily obtain the result that H 2 — ■> 1. 
From here we also obtain the result that if 7 > 4/3 that VLu — > 0. 



IV. QUALITATIVE ANALYSIS OF THE CASE = A = 

In the invariant set U + U U°, we can show that q > — 1. This implies that the invariant 
set f^A = is the invariant set containing all of the past asymptotic behaviour (i.e, early 
times) for all ever-expanding models in U + UU°. It can be argued that a scalar field becomes 
essentially massless as it evolves backwards in time, hence it will dominate the dynamics at 
early times (see Section ^). In our effort to understand the dynamical behaviour at early 
times we shall assume in the analysis that follows that there is no perfect fluid (Cl = 0) and 
the four-dimensional cosmological constant is zero (Cl^ = 0). 
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If Q = Q\ = then the equations (|3.3|) reduce to 



D' 


= -(l + q)HD 


(4.1a) 


H' 


= q(H 2 - 1) 


(4.1b) 


K 


= (q-2)HE+-S+ 


(4.1c) 




= (q - 2)#£_ - S- 


(4.1d) 


N[ 


= N x {Hq-A%) 


(4.1e) 


N' 


— AT ( fjn _i_ 9Y 1 j_ o./qv 1 "\ 


(A tf\ 


XT' 
N 3 


AT ( TJ ~ I or <~K /o"v \ 

= l\3{iiq + 21j + — zvdl- l -) 


//II «A 




a/6 ~ 


(4.1h) 


*' 


= (g-2)#$-^-M> 






(4.1i) 




= 2M u (q-l) 


(4-lj) 



The quantity g is the deceleration parameter, and S + and 5*_ are curvature terms that 
are now defined by the following expressions: 

q = 2t 2 + + 2t 2 _ + 2^ 2 - <§ 
'l 



5 4 



+a 

1 
6 



~D 2 (^r 2 + $)(5^ 2 -$) + Q W 



N 2 -N, 



S_ = ~ V3 (iV 2 - jV 3 ) (-iVi + iV 2 + iV, 



The two constraint equations are 
1 



H 2 + ^(N 1 N 2 N 3 ) 2 ^ 



h 2 = t 2 ,+t 2 _ + n total 



+ — [ (N 2 + N 2 + N 2 ) - 2(N 1 N 2 + N 2 N 3 + N t N 3j 



where 



n 



total 



.,2 Aotal 

3D 2 



$ + * 2 + a 2 



u 



(4.2) 

(4.3) 
(4.4) 

(4.5a) 
(4.5b) 

(4.6) 



A. Equilibrium Points at Finite D 



Here we define X = [D, H, £ + , Aq, A^, A^, f^] and we restrict the state space 
accordingly to be S = {X G = 0, = 0}. The equilibrium points X , can be classified 
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into one of the three matter invariant sets that do not have a perfect fluid component. 
Note, e is a discrete parameter where e = 1 corresponds to expanding models, while e = — 1 
corresponds to contracting models. The notation [xm] signifies that the preceding eigenvalue 
has multiplicity m. 



1. Vacuum, °il ' 

R e ; ROBERTSON- WALKER (Radiation) 

X = [0, e, 0, 0, 0, 0, 0, -y] The eigenvalues in the nine dimensional phase space [flu 
eliminated via Eq. ( |4.5b| )1 are 



e(-2, -1,-1, -1,1, 1,1,2,4). 

This point is obviously a saddle with a five dimensional unstable manifold (e = 1), 
and hence when the dynamics are restricted to the remaining constraint surface Eq. 



( |4.5a| )1, this point will remain a saddle (It can be shown that the stable manifold on 



the remaining constraint surface has dimension four) . 

R n e ; BIANCHI II (Radiation) 

X = [0, e, e|, 0, 0, 0, 0, 0, g^-] and two other equivalent points obtained through 
transformation (|3.8|) . The eigenvalues in the nine dimensional phase space [fly elimi- 
nated via Eq. ( [4.5b|) ] are 



e (-2,-l,-l,-i(l±v^,5 5 2,4) 

This point is obviously a saddle with a four dimensional unstable manifold (e = 1), 
and hence when the dynamics are restricted to the remaining constraint surface Eq. 
(|4.5a| )1, this point will remain a saddle. 



K e ; KASNER (Vacuum) 

X = [0, e, cos(6 l ), sin(6'), 0, 0, 0, 0, 0, 0] where — n < 9 < ir. This is a special case of the 
Kasner surface defined in next section. 

B VII e ; BIANCHI VII (Vacuum) 

X = [0, e, — e, 0, 0, s, s, 0, 0, 0] and two other equivalent lines of equilibria obtained 
through transformation ( |3.8| ) where 0<s<oo. Ass^O, these lines of equilibria 
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approach a point on K e . The eigenvalues in the nine dimensional phase space [flu 
eliminated via Eq. ( |4.5b )] are 



e(-3,0,0,2,4,6,6,±2si) 

One of the zero eigenvalues corresponds to the fact that this is a one-dimensional 
set of equilibrium points. The span of eigenvectors associated with the eigenvalues 
[0, 2, 4, 6, +2si, — 2si] does not include the '.D' direction. The eigenvector associated 
with the eigenvalue —3 is the only eigenvector with a component. Therefore this 
point is a saddle (e = 1) with a one dimensional stable manifold in the eight dimen- 
sional phase space S when the dynamics are restricted to the remaining constraint 
surface Eq. (|4.5a| ). 

2. Massless Scalar Field, °n°' ± 

/C e ; KASNER (Massless Scalar Field) 

X = [0, e, sin(<^) cos(#), sin(<£>) sin(6 l ), 0, 0, 0, cos(0), 0, 0] where — n < 9 < n and — | < 
if < |. Or + E 2 + \P 2 = 1. The eigenvalues in the nine dimensional phase space 
[Cly eliminated via Eq. ( [4.5b|) l are 

e(-3, 0, 0, 4, 4) and 2e - 4£ + , 2e + 2S + ± 2^3E_, 6e + V6k^ 

The two zero eigenvalues correspond to the fact that this is a two-dimensional set of 
equilibrium points. It can be shown (as in the B VII € case) that the eigenvalue —3 
corresponds to the D direction, therefore, this set of equilibria will always be a a set 
of saddle points with an unstable manifold having dimension no more than seven. 

3. Massive Scalar Field, °il + ' ± 

T+\ ROBERTSON- WALKER (Positive Curvature, Scalar Field) 

X = [0, ^fi, 0, 0, V4 - 2k 2 , V4 - 2k 2 , - 2k 2 , |, 0] The point only exists for 

k 2 < 2 and as k 2 — > 2 this point approaches P e . This is the only equilibrium point for 
finite values of D for which we can directly calculate the eigenvalues of the dynamical 
system restricted to the eight dimensional phase space S. The eigenvalues in the eight 



14 



dimensional phase space [H, Clu eliminated via Eqs. ( |4.5a| and |4.5b )] are 



e(-^fk, 0, -^(k ± v /F + 8 (P-2))[x2], -^(k ± Vfc 2 + 4(2-fc2))) 

This point when it exists (e = 1) has a six dimensional stable manifold, a one dimen- 
sional unstable manifold, and a one dimensional center manifold. 

P e ; ROBERTSON- WALKER (Zero Curvature, Power Law Inflation) 

X = [0, e, 0, 0, 0, 0, 0, — ^jp, 1 — 0] The eigenvalues in the nine dimensional phase 
space [Q u eliminated via Eq. (|4.5b|) 1 are 

e(~ k 2 - 2, \{k 2 - 2)[x3], \(k 2 - 6)[x3], k 2 - 4) 

If k 2 < 2 (e = 1), then this point is a local sink in the eight dimensional phase space 
S. If k 2 < 2 (e = — 1), then this point is a local source in the eight dimensional 
phase space S. When 2 < k 2 < 4 and when 4 < k 2 < 6 this point is a saddle in the 
nine dimensional phase space (e = 1), and hence when the dynamics are restricted to 



the remaining constraint surface Eq. ( 4.5a ), this point will remain a saddle. When 



k 2 = 4 this point experiences a bifurcation with the point TZSJ 7 ®, the Robertson- 
Walker radiation-scalar field scaling models. When k 2 = 6 this point becomes part of 
the Kasner massless scalar field models K e . 

RSF° e ; ROBERTSON- WALKER (Zero Curvature - Radiation-Scalar Field Scaling Model) 
X = [0, e, 0, 0, 0, 0, 0, 3J2, pjjr] When k 2 > 4 this point is an element of U + and 

when k 2 < 4 this point is an element of U~ . The eigenvalues in the nine dimensional 
phase space [flu eliminated via Eq. (|4.5b[) 1 are 



e(-2, -1, -1, 1, 1, 1, 2, -±-(k ± ^ + 16(4-^))) 

This point (e = 1) is easily seen to be a saddle with a four dimensional unstable 
manifold when k 2 > 4, and a five dimensional unstable manifold when k 2 < 4, and 
hence when the dynamics are restricted to the remaining constraint surface Eq. ( |4.5a|) , 
this point will remain a saddle in all cases. Note, the cosmological model represented 
by this equilibrium point has the property that the energy density attributed to the 
scalar field is proportional to the energy density of the dark radiation coming from the 
bulk 



15 



RSF n e ; BIANCHI II (fl u + 0), Radiation-Scalar Field Scaling Model 

X = [0, e, |, 0, 0, 0, — gp-, 7 gfc2~2 2 ] and two other equivalent points obtained 
through transformation ( |3.8| ). The eigenvalues in the nine dimensional phase space 
[flu eliminated via Eq. ( [4.5b| )1 are 



-2-1**2 

' '2'2' ' 



2k 2 ± yj (2k 2 ) 2 - 2k 2 (23k 2 - 64) ± y/(23k 2 - 64) 2 - 64k 2 (7k 2 - 32) 

This point is always a saddle point in the nine dimensional phase space, and hence 
when the dynamics are restricted to the remaining constraint surface Eq. ( [4.5a]) , this 
point will remain a saddle in all cases. Note, the cosmological model represented by 
this equilibrium point has the property that the energy density attributed to the scalar 
field is proportional to the energy density of the dark radiation coming from the bulk 

p sf = ^((p) 2 + V(0)(xU = pu. 

SF n e ; BIANCHI II (fl u = 0), Scalar Field Model 

Xo = [0, e, =^P, 0, ey ^ggg , 0, 0, -g$g, 36^1^, 0] and two other equiv- 
alent points obtained through transformation ( |3.8| ). This point only exists when 
2 < k 2 < 8. As k 2 — > 2 this point approaches the point P e and as k 2 — > 8 this 
point approaches a point in the set /C e . The eigenvalues in the nine dimensional phase 
space [flu eliminated via Eq. (|4.5b|) 1 are 

6(k 2 - 8)[x2], 16(A; 2 - 2), 12(k 2 - 2)[x2], 2(7k 2 - 32), -9k 2 , 



k 2 + 16 

3[(k 2 - 8) ± v / (A; 2 -8) 2 + 12(A; 2 -2)(A; 2 -8)] 

This point is always a saddle point in the nine dimensional phase space, and hence 
when the dynamics are restricted to the remaining constraint surface Eq. ([4.5a]) , this 
point will remain a saddle in all cases. The cosmological model represented by this 
equilibrium point is that of a Bianchi II model with a exponential potential scalar 
field. 
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Note there are other equilibrium points to the dynamical system ( |4.1| ), but they corre- 
spond to points that are neither inside the set S nor on its boundary, and do not represent 
any limiting behaviour of the Bianchi Type IX brane-world scalar-field models. 

B. Initial singularity 

In the analysis of local equilibrium points at finite D we found in the previous subsection 
that there were no local sources corresponding to expanding models. In order to complete 
the analysis, and determine the dynamical behaviour close to the initial singularity, we need 
to examine what happens as D — > oo. In this subsection we present a heuristic analysis of 
the situation, and include ordinary matter, that is O 7^ 0. 

We define a new bounded variable 

d =p^l> -l<d<l (4.7) 



and examine what happens as d — > 1 (assuming H > 0). From Eqs. ( |3.3a| and |3.3b| ) (or 



( [4. lap and |4.1b| )) we have that d — > 1 and H — > 1 monotonically (see also Section |V|) and 



hence we need to consider the equilibrium points in the set d — 1. 
The analysis depends on whether the quantity 

,2 



A= ^(fi + ^ 2 + $) ((3 7 - l)f2 + 5^ 2 - (4.8) 



in the expression for q in equation fl3.4|) is zero or not in an open neighborhood of the 
singularity. If A 7^ 0, and assuming A > (see Section 0), we define a new time variable 
by 

AH 

and the remaining evolution eqns (on d — 1) become 

S+ = E+, S_ = S_, N^Nx, N 2 = N 2 , N 3 = N 3 , 

n = 2C1, n A = 2n A , ^ = i = 2$, n u = 2Vt u . (4.9) 

Therefore, the only equilibrium point is 

£+ = S_ = N x = N 2 = N 3 = n = n A = * = i = VL U = 
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and this is a local source. This equilibrium point corresponds to the Brane- Robertson- Walker 
solution [O with 

— 2 

^-D 2 ^ 2 = 1, _D 2( J = (D 2 f2 = 0) 

(see Section |V|) . 

It remains to consider the case A = 0. However in this case 

q = 2(t 2 + + t 2 _) + 2^ 2 - $ + 37 ~ 2 Q - fi A + a 2 fi w (4.10) 

and resulting dynamical system ( |4.1b|) - ( fl.lj ) has the same local equilibrium points as in 
the previous subsection but with d = 1 (instead of d = 0). In particular, there are no local 
sources, and the Kasner equilibrium points are saddles [0. Hence all orbits asymptote in 
the past to the BRW source. This will be discussed in more detail in the next section. 

C. Local stability of the Brane Robertson- Walker solution 

In order to analyze the dynamical system ( f4.1|) for large values of D, we can define the 
following new variables 

$ = r 2 sin 2 9, and ^ = r cos 9, Q x = ^a 2 D 2 {$ + ^> 2 ) 2 = \ 2 D 2 r A 

Using these new cylindrical coordinates, the infinite variable D is essentially replaced by 
the bounded variable Q\ in the set U UU°, hence in this set the only variables that remain 
unbounded are the NJs. 
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The dynamical system ( |4.1|) becomes 

H' = q(H 2 - 1) (4.11a) 

S' + = (q-2)Ht+- S + (4.11b) 

E'_ = (q- 2)#E_ - 5_ (4.11c) 

N[ = 7V!(#g-4E + ) (4. lid) 

N' 2 = N 2 (Hq + 2t + + 2V3t_) (4.11e) 

JV£ = N 3 (Hq + 2t + -2^3E_) (4.11f) 

^ = 2#Q w (g-l) (4.11g) 

ft' A = 2im A (g + l-6cos 2 #) (4.11h) 

r' = rH(q - 2 + 3 sin 2 6) (4.111) 

/fi 

0' = sin0(3#cos0 + — Ax) (4.11j) 

where S + and >!?_ are curvature terms defined previously and where 

q = 2E* + 2E 2 + r 2 (3 cos 2 - 1) + fi A (6 cos 2 - 1) + a 2 fi w (4.12) 

The two constraint equations become 

1 = H 2 + ^(N 1 N 2 N 3 ) 2 / 3 (4.13a) 
H 2 = t 2 + + t 2 _+Q x + a 2 (l u + r 2 

+ ^ ((iV 2 + iV 2 + N 2 ) - 2(iV x iV 2 + N 2 N 3 + iViiVg)) (4.13b) 



From the constraint equations, we have < H 2 , T? + , E 2 ,, flu, Q\, r < 1, < < 7r, and 
NiN 2 N 3 < 8 (assuming ft w > 0). 

The Brane-Robertson Walker solution |I(J is represented by an equilibrium point in the 



set tt x ^ (D -> oo). If we let X = [H, E+, E_, JVi, N 2 , JV 3 , fl u , fi A , r, 0] then this new 
equilibrium point is 

m e ; Brane Robertson Walker 

X = [e, 0,0, 0,0, 0,0,l,0,f ± §] and X = [e, 0, 0, 0, 0, 0, 0, 1, 0, f ]. Using the con- 
straint equation to eliminate fly, the eigenvalues of the linearization at the point 



9 = 0, 7r are 



e(10,10,5,5,5,3,3,3,3) 
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and when 9 = | the eigenvalues of the linearization are 

£(-3,-3,-3,-2,-2,-1,-1,-1,0) 

The values 9 = 0, ir correspond to invariant directions in the neighborhood of this equilib- 
rium point. For example 9 = 0, tt corresponds to $ = invariant surface (that is the massless 
scalar field models). The value 9 = n/2 corresponds to \£ = is not an invariant direction, 
and we note that 9' > in a neighborhood of the equilibrium point near \l/ = 0. We easily 
observe that this equilibrium point is a source that strongly repels away from \1/ = 0. That 
is, if one was to approach this equilibrium point along a typical orbit (backwards through 
time), one would asymptotically approach a massless scalar field Brane Robertson- Walker 
solution. 



TABLE I: Local sinks and sources. Note that each of these local sinks has a corresponding local 
source, P- and m_|_. 



Sink 


Conditions on k 


Description 


P+(Zero curvature, power law inflation) 


k < \[2 


D = 0, H = 1, expanding models. 


m_ (Zero curvature, isotropic Braneworld) 


all k 


D — > oo, H = —1, contracting models 



V. GENERAL DYNAMICS 

We now return to the general case (i.e., equations ( |3.3| ) - ( |3.6| )) and include a general 
potential for the scalar field, and study what happens as D — > oo at the initial singularity 
(assuming this occurs). We shall also assume normal matter with 1 < 7 < 2. 

We recall that there are no sources for finite values of D. Also, for a general scalar field, 
equations (|3.3jD and (|3.3k| ) become 



^' = (q-2)M-e§ (5.1) 
<§' = 2(q + 1)M + 2e$^ (5.2) 

where e (related to the usual inflationary slow roll parameter e) is defined 



"2kV 
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From the Friedmann equation we have that 

n 2 

-D 2 (Q + $ + ^ 2 ) 2 = Q x < 1 (5.3) 



4 

and hence each term DQ, D§, D^> 2 is bounded (since the left-hand side is the sum positive 
definite terms). Hence, as D — > oo, Cl, \I/ 2 , $ — ► 0. 

It is easy to show that Q u —>■ as D — ► oo. Hence fl3.4| ) becomes 

2 

g = 2(£ 2 + £ 2 ) + ^-L> 2 ([fi + ^ 2 + 4>][(3 7 - l)tt + 5§ 2 - $]J (5.4) 
Assuming H > 0, equations ( |3.3a| ) and (|3.3b|) imply that as r — ► — oo for D — > oo, either 



q — > or q is positive in a neighborhood of the singularity (g can oscillate around zero; 
indeed, it is the possible oscillatory nature of the variables — i.e., q need not be of single 
sign — that causes potential problems). However, if q — > 0, Q3.3h|) implies 



(V = HVt{2 - 37) (5.5) 

which implies a contradiction for 7 > | (i.e. Q -/+ as r — > —00). Hence, as r — > —00, 
q > 0, where D —>■ 00 and 

= _qf(i - # 2 ) (5.6) 

and hence if — > 1 (assuming positive expansion) monotonically. [Note: this implies the exis- 
tence of a monotonic function, and hence there are no periodic orbits close to the singularity 
near the set H — 1 — all orbits approach H = 1.] In addition, (ft.7a ) gives (iV 1.^2^3) ~~ ^ 0. 
From Eqs. ( |3.3h|) , ( |3.3i| ), ( |3.3j| ) we have as r — > — 00 

■^--►0, (5.7) 

that is, the scalar field becomes effectively massless, and 

£- if 7 < 2. (5.8) 

This follows directly from the evolution equation in the case of an exponential scalar field 
potential, e = \/%k, and follows for any physical potential for which e is bounded as r — > 



g = 2(S^ + S 2 ) + (3r- 1)C 2 (5.9) 



2 

-00. Hence we obtain 
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where T = 2, C 2 = ^D 2 ^ 4 if 7 < 2; T = 2, C 2 = ^D 2 {Q + ^ 2 ) 2 if 7 = 2; T = 7 and 
C 2 = \D 2 n 2 if there is no scalar field, where VL\ — > C as r — ► —00, and the Friedmann 
equation becomes 

1 - {t 2 + + £ 2 ) - C 2 = -^((iV 2 + N 2 + N 2 ) - 2(^x^2 + N X N 3 + N 2 N 3 )) > (5.10) 



From Eqs. ( p.3h|) , (|3.3j|) we obtain 

(C 2 )' = 2C 2 H[q-(3T -1)}. (5.11) 

We still have the possibility of f}± or oscillating as the singularity is approached (as in 
the Mixmaster models). A rigorous proof that oscillatory behaviour does not occur can be 
presented using the techniques of Rendall and Ringstrom |27| and |]28| (using analytic ap- 



proximations to the Brane-Einstein equations for E^, N a ; i.e., estimates for these quantities 
that hold uniformly in an open neighborhood of the initial singularity). We can then prove 
that S-t — >• as r — > —00, Cl\ — > 1, and we obtain the Brane Robert son- Walker source. 
Alternatively, from equation ( 5.11|) we have that either C 2 — > or q — > 3T — 1. If C = 



then g = 2(E^ + S 2 ) and we can show that we obtain contradiction (C = implies D = 0!). 
Hence C 2 7^ 0, and g — * (3r — 1), so that 

t 2 + + t 2 _=$ (5.12) 

and it is a simple matter to show that iV Q — > and again we obtain the Brane-Robertson 
Walker source. 

It was shown earlier (for no perfect fluid and a scalar field with an exponential potential) 
that the Brane-Robertson- Walker solution is always an equilibrium point of the system and 
that local stability analysis shows that it is a local source. 



VI. CONCLUDING REMARKS 

Assuming that Vtu > 0, the future asymptotic behaviour of the Bianchi IX brane world 
containing a scalar field having an exponential potential is not significantly different than 
what is found in general relativity p9| . We observe that for < k < y/2, the future 
asymptotic state is characterized by the power-law inflationary solution, and if k > \pl 
there no longer exists any equilibrium point representing an expanding model that is stable 
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to the future. We therefore conclude that if k > a/2 then the Bianchi IX models must 
recollapse. In |2{J it was shown that if k > v2 then a collapsing massless scalar solution 
was a stable equilibrium point. Here, in the Braneworld scenario, we have that this final 
end-point is the Brane-Robertson Walker solution. However, we also observe that a typical 
model on its way to this final endpoint will asymptote towards a collapsing massless scalar 
field solution. However, if Cly < 0, then a variety of new behaviours are possible including 



possible oscillating cosmologies [|i~2 



The past asymptotic behaviour of the Bianchi IX braneworld containing a scalar field (and 
ordinary matter) is significantly different than what is found in general relativity. It is known 
that the Bianchi IX perfect fluid models approach a MixMaster attractor (Kasner saddles 
joined by Taub separatrices) towards the past and are known to have chaotic behaviour. 
Here we observe that the Brane Robertson- Walker solution is a global source, and that 
there is no periodic behaviour near the initial singularity. 

Due to the quadratic nature of the brane corrections to the energy momentum tensor, 
a rich variety of intermediate behaviour is possible in these two fluid models. We note the 
existence of the dark radiation density U together with a scalar field is similar to previous 



analysis done in |30[ on scaling models. Here the equivalent equation of state is pu = (§ — 

It is known that the bifurcation value for these scaling models is k 2 = 37, which for 7 = 4/3 

corresponds to a value k 2 = 4, a bifurcation value found in the analysis above. 

The intermediate behaviour of these multifluid models considered here is extremely com- 
plex, and is not discussed in detail here. However, a more complete analysis (in which some 
of the intermediate behaviour is outlined) of the Bianchi type II models is currently under 



investigation 31, 32 
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Appendix: Integrability conditions for V^u = 0, Q M = 

The complete set of four dimensional Brane equations for a general imperfect fluid energy 
momentum tensor is given in Maartens (see equations 26-29 and A1-A10 in Q). If = 
0, <2 M = 0, then in the case of a perfect fluid with an equation of state p = (7 — l)p, we 
obtain the following equations when the vorticity is zero (u = 0): 



(7-l)l n 



7 i> 

p + 7 9p = , 

w + lew = , 

3 

6 + ie 2 + - D"^ + + ^ 2 (3 7 - 2)p - A 



6 



U 



DV 



A^A v) 



„„ -!)„(-) 



-K 



V» 3 --^r L -,- JH 3A 

and the Gauss-Codazzi equations on the brane 



K 2 \ 



1 

l (iiv) + 3 < " 
2 



ir + -B 2 - 0^0^ - 2k z P - 2A 
3 

Using Eq.( |6.1aj) , Eq. (|6.1d|) becomes for 7 7^ 1: 

"4(7-1)1/ K «7 



! „ 12 



P 



IV- 



(6.1a) 
(6.1b) 
(6.1c) 

(6.1d) 

(6.1e) 
(6.1f) 
(6.1g) 
(6.1h) 



(6.1i) 
(6-lj) 



(6.2) 



7 p 2 

Taking the directional time derivative (i.e., a dot) of Eq. (|6.2|) , and using Eqs. ( |6.1c|) and 
( |6.1d|) , and using known relations for interchanging space and time derivatives, e.g., 



VPJ]' = h^lVu + A v ]f - 



(6.3) 
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we obtain the integrability condition 

^(37-4)W-^-p 2 
4 



pD„0 



4 

k 7 



— (3 7 - 4) ( 7 - 1)W + ^(37 - l)p 2 
37 6 



OD^p. 



(6.4) 



Taking the directional time derivative of ( |6.1g| ), interchanging space and time derivatives, 
and using Eqs. fl6H) , (|6Uj ), (l6lh|) (and (|6Ub|) , flOc?) and ( grp ), we obtain 



^ ^[12W + K 4 p 2 ]iD M p = 0. 

7 p 

Hence for 7 7^ 1 and [12W + /t 4 p 2 ] 7^ 0, we must have 



(6.5) 



D M p = 0, 



so that 



A M = T)fjl4 = D M 9 = DV p = D^o" = 



and hence in general the brane is spatially homogeneous. 
The special case 

7 = | 12W + K 4 p 2 = 0, 



(6.6) 



in which the integrability conditions yield no constraints, corresponds to the case in which 
there are no corrections to the general relativistic equations. 
Alternatively, Eqs. ( |6.1b| ) and (|6.1c| ) imply that 



U = hp^; h = 



(6.7) 



and so Eq. ( |6.2j ) becomes 



Defining 



and a 



2(3 7 -4) 
(3 7 -2) 



and b 



3«V 
4(37-2) ' 



1 2(3 7 -2) 
\$) = -n 3 7 

h 

we then obtain (for 7 7^ 



which integrates to 



# (l-o)[(l-a) 



(6.8) 



(6.9) 



(6.10) 
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where D u f = and a = -r^— 



4^4-37) • ^his 10| ) essentially yields the functional 

dependence for p. In the exceptional case 7 = 4/3, p is separable in local coordinates. (From 
Eq. (|6.8|) it can be seen that the case 7 = 2/3 is also exceptional). 
Eq. ( |6.4| ) then becomes 

-4, 



'4 ftS 2 n 4 

1(37-4)/.-^-^ 



4 (3 7 _4)( T _1) K 4 7(7 _ 1/3) 2 

n p 4 t 

37 2 y 



P 

(6.11) 

From Eqs. ( |6 . 1 0|) and (|6.1cj) , and using local coordinates so that h = /i(x 7 ), / = f(t), 
we obtain the functional forms for p and G, whence on substitution into Eq. (|6.11| ) (for 
7 7^ 4/3) we obtain a contradiction if D M p 7^ 0. 

Finally, in the special case 7 = 1, in which p = and = 0, we obtain analogues of Eqs 
3) and (|6lp (which do not follow for = 0): 

1 



D M W = --« pD^p 



e D « e 



-K 4 pD M p. 



(6.12) 
(6.13) 



In local coordinates we then find (for D^p 7^ 0) 



1 + ^4 2/3 + 1,8 4/3 

8 H 32 H 



(6.14) 



which implies that p = (and hence = 0)! 
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